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- summary of results 


his article summarizes laboratory results from exper: 

iments conducted on a rotor model. These experi- 

ments proved that a model of Fluid dynamic forces 
acting on a concentric, or slightly eccentric, clastically-sup- 
ported cylinder rotating in a stationary cylinder can ade- 
quatelv predict the system natural frequencies of the first 
and higher modes, instability thresholds and the existence of 
self-excited vibrations. Seleexcited vibrations. which are also 
known as whirl or whip, occurred, due to interaction with 
the surrounding fluid. This article also discusses the physical 
factors that affect fluid whirl and whip. 


Introduction 

The dynamic characteristics of most fluid-handling 
machine rotors depend on physical phenomena taking place 
in rotor-surrounding fluids. In particular, the phenomena 
occur in interstage seals, blade-tip clearances and fluid-buhri- 
cated bearings, In all these cases, a description of the physi- 
cal situation can be simplified as a cylinder rotating inside a 
stationary cylinder with a relatively small clearance between 
them Independently from the axial component of the flaw 
inside the clearance, the rotating cylinder (rotor) forces the 
fluid, through viscous friction, to rotate, creating a circum. 
ferential flow. The rotating fluid begins to participate in the 
system dynamics. New forces are generated in the fluid film, 
which, in a feedback loop, act on the rotor. The rotor 
responds with lateral vibrations, resulting in a forward, 
almost circular precession, known as “Пи whirl" апа 
"whip" [1]. 

The deseription of £luid whirl and whip that follows is 
based on the example of a flexible rotor supported by a lat- 
crally rigid pivoting brass oilite bearing at the inboard end, 
and by a 360 oil-lubricated bearing at the outboard end. In 
order to set the rotor journal at any radial position inside the 
outboard bearing clearance, adjustable radial supporting 
springs were attached to the rotor through a rolling clement 
bearing (Figure 1). The concentric journal position was cho- 
sen as the first step, because it creates the most favorable 
conditions for a fluik-induced instability to occur. 


Fluid whirl and whip 
Rotor lateral vibrations are measured by two sers of XY 
proximity probes (Figure 2). Data in the upper plot is from 


Orbit £ 


probes mounted at rotor midspan; data in the lower plot is 
from probes mounted at the Fluic-lubricated bearing, Each 
full spectrum cascade plot includes Plus Orbit plots, which 
show the magnified actual paths of the rotor centerline at 
various rotative speeds. At low rotative speed, the only later- 
4l vibration component in the spectrum is svnchronous 1X 
vibration, due to an unbalance in the rotor (Figure 2). The 
rotor unbalance acts as a centrifugal rotating force which 
excites rotor 1X responses. A mew vibration component 
appears in the spectrum at the rotative speed known as the 
instability threshold (onset). n machines, such as the onc 
considered, it occurs at a speed that is usually lower than the 
machine's nameplate first balance resonance (critical At the 
threshold of instability, the rotor motion becomes unstable, 
and the lateral vibration amplitudes start increasing. At First, 
this growth is almost instantaneous. Then it slows down, duc 
to the nonlinear effects in the fluid film at higher eccen- 
tricity. In. a few rotations, a new state of equilibrium is 
achieved. This new component in the rotor vibrational 
response spectrum is the fan! cycle of self-excitee vibre- 
ows, known as fifa woe or in this case, aff aod The fre- 
queney of fluid whirl is lower than rotative speed, £2 , and, 
for the considered system is 0,475 £2 {їп general, it can be 
апу fraction, A). The important characteristic of the fluid 
Whirl is that its frequency is proportional ta the retative 
speed. This proportionality can be observed on the speg- 
trum cascade plots (Figure. 2) between 2400 to 26000 rpm 
and 4300 to 6000 rpm. 

During fluid whirl, che mode of rotor vibration is conical, 
like a rigid body, The journal orbits are circular, and forward. 
When the rotative speed approaches the machine's first bal- 
ance resonance, (the rotor lateral moide first natural frg- 
queney), the LA vibration amplitudes start increasing. The 
rotating unbalance force pushes the journal to the side or 
the bearing, resulting in high dynamic cecentricitv. This 
action causes the flow pattern to change and the fluid film 


racial stiffness to increase, This action results in a temporary 
suppression Of the fluid whirl [2], changing the flow pal- 
tern. and inereasing the fluid film radial stiffness [2]. The ces- 
sation of the rid whirl follows, In this case, the unbalance 
farce acts similar to a radial constant force (such as gravity, 
tor instance, which causes the journal to move to a higher 
eccentricity, This modification of the flow pattern, anc 
increase of the radial stiffness, stabilizes the rotor The fluid 
whirl ceases, In contrast to a unilateral radial force, the 
unbalance force acts in a rotational fashion, not in a con- 
stant direction and results in 1X dynamic eccentricity, not a 
static one, The cfteets are, however, similar. The rotative 
рес width af the rotative speed range, where there is no 
whirl, depends on the amount of rotor unbalance, Por a bal- 
anced rotor, it is reduced te zero, 

lt i$ important to note that the mechanical resonance Ire- 
quency is not Constant it changes as à function of the jour- 
nal eccentricity ratio, At à low eccentricity гапо, where the 
journal is located near the center of the clearance, the fluid 
him radial stiffness and mechanical resonance frequency are 
lower At à high eccentricity ratio, where the shalt is posi- 
tioned near the periphery of the clearance, the Fluid film 
radial stiffness and, therefore, the mechanical resonance Ire- 
queney are higher. As a matter of fact. this resonance fre- 
queney corresponds to the case of a rigidh-supported rotor. 


The high "eccentricity" resonance frequency is typically the 


manufacturers “nameplate critical." The fluid-incuced reso- 
nance frequeney is also not constant: it changes as а function 
af the eccentricity ratio. At a high eccentricity ratio, the fluid 
average circumferential velocity ratio A is low and the fluid- 
induced resonance frequency is lower At a low eccentricity 
ratio, the fluid average circumferential velocity ratio CA) is 
high and the fluid-inducced resonance frequeney is higher, 
When continuing the rotative speed increase, after the 
first balance resonance, 1X amplitudes decrease, the forced 
dvnamic eccentricity of the journal is reduced, flow pattern 





Figure 1 


Rotor/bearing system. A - speed controller, B - 75 W electric 
motor, € - speed controller transducer, D = elastic coupling, 
E - Кеурһаѕо (once per turn) transducer, F - laterally 
rigid, pivoting bearing, G - 9.5 mm (0.375 in) diameter steel 


shaft, H, K - four proximity eddy current transducers, 
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mounted in XY configuration, I- rotor disk of mass 1.63 kg 
(3.616 Ib), J - four radial supporting spring systems, L- oil 
(T-10) lubricated bearing with 51 mm (2.0 in) length, 22 
ym (8.6 mils) radial clearance and two-port supply oil with 
10342 Pa (1500 Ib/in?) pressure, 
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Figure 2 


Full spectrum cascade of the rotor midspan (a) and journal. 


(b) lateral vibrations with orbits at different rotative speeds. 
Using information from two orthogonal transducers, the full 
spectrum displays forward and reverse circular components 
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of rotor elliptical orbits of all frequency components, The 
orbits at $370 rpm show the whirl pattern. Orbits at 6770 and 
8070 rpm show the whip pattern. 
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goes back to a predominant circumferential Mow, and the 
fluid whirl reoccurs at the second onset of instabili The 
whirl characteristics are similar to those in the lower rota- 
tive speed range, maintaining the same proportionality fac- 
tor of its frequency to the rotative speed. When the rotative 
speed approaches a value that is approximately twice the 
first balance resonance, the whirl frequency no longer main- 
tains proportionality to the rotative speed. However, the fre- 
quency ratio starts decreasing. (In the considered case, this 
rotative speed is actually the speed close to the inverse of 
(A) times the first natural frequency, 6, , 1/0475 = 2.1104, 2. 


While the frequency ratio decreases, the whirl ampli- 
tudes. especially at the rotor midspan location, signilicantly 
increase, With the rotative speed continuing lo Increase, al 
the end of this smooth transition, the frequency of the sell- 
excited) vibrations asymptotically approaches the constant 
Frequency which corresponds to the rotor high eccentricity 
first bending mode natural frequency, ey, . This rotor self- 
excited beibrgdion is called Fluid whip cos ер in this case. 
The rotor mode is no longer a conical, rigid body. It is the 
rotor first bending mode, which is associated with i . The 
whip orbits at the journal are circular, borward, and their 
amplitudes are almost. equal to the bearing clearance. At 
midspan, they are much larger and slighily elliptical, since 
the rotor system has some residual lateral anisotropy. The 
unfiltered orbits in Figure 2 slow the effect of two major 
frequency component interactions: The self-excited and the 
ІХ vibrations. The small black «ots in the orbits plots are 
generated by the Kevphasor transducer, and they indicate 
the vibration-to-rotation. frequency relationship. Two con- 
secutive dors indicate one shaft rotation. Therefore, the whirl 
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orbits of frequency 0475X will show two dots with the 
angle between them equal to 0.475 360° = 1717. COn the 
oscilloscope screen set in orbital mode, these two dots will 
appear to rotate in che direction against rotation. as at cach 
two rotations, they lag the orbit full circle by 3607 - 2 x 1717 
=18"), 

Fluid whirl and whip are two variations of the same phe- 
namenon, which is known as the fluid-induced instability of 
rotors. Two different terms arc used because the same phe- 
nomenon was discovered independently on two different 
machines: “whirl” in machines with an operating. speed 
lower than twice the first balance resonance and “whip” in 
high speed machines with an instability threshold. karger 
than twice first balance resonance. Since whirl and whip 
were observed independently, nobody, at that time, connect- 
ed them, 


A new natural frequency of the rotor/bearing system 
Fluid whirl and whip are both seltexcited vibrations. It is 
well-known that the selfexcited vibrations, as limit cycles of 
transition from an unstable mode, are conditioned by system 
nonlinearities. The limit cvele frequencies. аге, however. 
weakly sensitive to nonlinearities, and are equal to the nat- 
ural frequency of the system at the instability onset. There is 
no doubt about such a relationship in the case of the Fluid 
whip its Irequency corresponds ro che rotor high eccentric 
ity natural frequency of the first bending mode or nameplate 
critical. The following conclusion for fluid whirl is obvious: 
The whirl frequency, AG, must be close to the natural Fre- 
quency of the тогот bearing fluid system. In contrast to the 
classical natural frequencies of mechanical systems, which 
are based upon stiffness and mass, the whirl natural Fre- 
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Figure 3 


а) An example of fluid circumferential velocity ratio.A, ver- 
sus rotor eccentricity ratio esre (c ^ radial clearance) for 
two values of rotative speeds, obtained from impact testing 
of an operating rotor/bearing svstem [7]. (b) Spectrum ana- 
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lvzer screen display during impact testing. Three peaks 
appear: AQ, rotor bending mode natural frequency, and the 
rotative speed frequency, £2, 
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queney is ETuid-related, and its dominant characteristic is ris- 
cous damping. The radial damping in the fluid film has an 
over-critical value. There is also an important difference in 
the fluid film damping farce, in comparison to classical dash- 
pot camping. In the case of seals and lighth-loaded bearings, 
the fluid damping force rotates. The angular velocity of the 
fluid damping force rotation is AG. This rotation generates 
an additional, angular velocity-related damping force, which 
acts on rhe rotor in the fevieentivd direction. This is the dri- 
ving agent For the fuid whirl and whip. 

The analvsis of rotor/fluid svstems has generated the fol- 
lowing important results: 


(a The extension of Modal Analysis to include Cluid/solid 
interaction effects [3-6] 

(bo The adjustment of the fluid force models in seals andl 
lighth-loadedd bearings [7,8] 

(cO The extension of the Пиш! force model for other cases, 
such as blade tip clearance and rotating stall cases [9], 


Model of fluid film forces in seals and lightly-loaded 
bearings 

The fluid force model discussed below сап also be 
applied to seals and lighth-loaded bearings when an almost 
concentric cylinder rotates inside à stationary cylinder, with 
a small clearance, and when the dominant fluid flow pattern 
is rotor rotation and triction-induced circumferential flow. 
Por the concentric "cylinders? the average dynamic force of 
the Huid film in the radial direction, Г, expressed in the coor- 
dinate system xaJ, whieh rotates at the angular velocity 
AU, is as follows [7,8]; 


-Е = М,ғ, + Dr, + |K, tK (|r, |, e 
r =x, + fp, j= far [= dx “+ y 








where My, is Fluid inertia effect, D is fluid film radial damp- 
ing, and Ky and Rp are linear and nonlinear parts of the fluid 
film radial stiffness, respectively. The nonlinearity here is 


considered in the general form of a continuous function of 


the rotor lateral displacement г.) limited by the radial clear- 
ance, c In general, the fluid film stiffness is an increasing 
function of eccentricity: Eq. 1 yields the formal definition of 
AQ whieh can be found in Eg. CD it is an шийит velocity 
at weh Hoe fluid force rotates, Based on circumlerential 
flaw pattern analysis [10] A is called fuel circumferential 
average velocity ratio, The transformation тоге of Eq. 
(1) inte the stationary coordinates x, р vields: 


-Fz-EeU = М,(#—2 ДО» - XO'r]4 (2) 


+ Dr = ЎА 4 |K je ( zin = X Јр 


where F is the fluid force in stationary coordinates. In this 
iarmat, the Fluid force can be directly introduced into the 
classical rotor model. This force is acting on the rotor. In Ee. 
(2), x and y are lateral orthogonal displacements of the 
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гог (at the fluid force axial location) described in the sta- 
Honarv coordinate system. From the rotor stability stand- 
point, the most important component of the fluid force (2) 
is the tangential component, (OA (The "7" indicates that 
this term is perpendicular to the radial direction, à 3. The 
coefficient DAR is known as the “cross stiffness" term. In the 
model (25, chis "coefficient" is obtained based on flow aver- 
age velocity consideration, and with predominantly circum- 
ferential Mow, it results from the damping force rotation. It 
is. therefore, a producer of well-defined physical parameters 
D, A and ©. 

Ihe fluid force model (23 has been identified through 
extensive modal testing [4-6]. The simplified Пиш force 
model in which А = 0.5 has been mentioned in engineering 
literature since the early 1960s [11, 12]. However, this is only 
а тегу rough approximation. With an assumed constant 
value of A, the model cannot reflect important features of 
the fluid force. The ratio, A, is actually not a constant coeffi- 
cient, but a function of rotor eccentricity (Figure 4), rotative 
speed, and a function of the fluid parameters and Now pat- 
tern. Fhe Tow pattern can casily be modified by bearing or 
seal Пи inlet/outlet conditions, such as preswirls, geometry 
of operating “cylinder” surfaces (c.g, multidobe bearings. 
labyrinth seals), swirl brakes, fluid antiswirl injections 
(Figure 4 [13, 14, ete. The modifications of the Huid flow 
pattern, from a pure concentric rotor rotation-induced cir- 
cumferential pattern. inte a predominantly: axial pattern, 
reduces the value of A and improves the rotor stability. In the 
next section, it will be shown how À is related to the stabil- 
itv of rotor/bearing/seal systems. 


ROTOR 


PROCESS 
FLUID 
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Antiswirl injection principle: fluid injected tangentially to 
ihe seal in the direction opposite to rotation produces fluid 
torque which reduces rotation-cenerated fluid circumfer- 
nenial average velocity, AL? [14, 14]. 


Rotor/seal rotor/bearing and models 

The model of the rotor system with the assumption on 
system isotropy, limiting considerations to the first lateral 
mode and using Eg. (2) with fluid inertia neglected for clar- 
Iv}, is as follows (Figure 5x 
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Mr, + Dy HIK + K ,]r - K,r, = 0, (3) 
(е) xt) y(t). $212 


Dr. - jar, ) т |. +A, +A, + K „|, Jr. — K,n-ü 


where MD. are rotor modal mass and damping. respective- 
ly. Ky Kok, аге rotor partial modal stiffnesses for the 
rotor/scal model (Figure. бау Ky is the rotor. supporting 
spring stiffness for the rotor/bearing model (Fig. 5b. гү and 
r; are rotor lateral displacements at the disk and at the fluid 
foree location in the complex number format. 

First the svstem cigenvalucs will be discussed, The char- 
acteristic equation for the lincar system (5) (Ky, 200. is as fol- 
lows: 


[K, +K, +50, + Ms! | Ds - JA) + 


» (4) 
+K,+K,+K,|—K; =0 
where s is the eigenvalue. Introducing the complex eigen- 
value б) by s = fio, Eq. CH can be transformed inte the fol- 
lowing approximate Гот: 
| D.K; | 

DK, +K, -MXQ Y + KES] 
| 2 | 





H к -MY | 
+2 K, +X, T E 
Di" * ^CK--K,-MXQ' | 


where on the right side of Eq. (5) 0 was replaced bv А52, (the 
zeroth approximation), and some terms with D, were omit- 
ted. The real part of Eq. (55 represents the first approxima- 
tion to one of the system natural frequencies. Since the 
second term of the real part of (5) is proportional to the 
rotor damping, D., its value is small. Therefore, this frequen- 
cy is approximately AQ. The imaginary part of Eq. (5), when 
equalized to zero, provides the instability threshold £2: 


о„=1 къ EK) т 
i A \ M | K 2 + (K, " К, | | 


(a) 


nt) rq (f) 


The last expression under the radical in Eq. (6) represents 
{мо stilfnesses, Ky and КК in series, Therefore, it results 
in а stiffness lower than either of the components results, In 
the rotor/seal case {Figure 5a), this produces a value which 
may significantly contribute to Ку. As a result, the instability 
threshold (63 becomes relatively high. often as high as twice 
the first natural frequency of the system. In this range, the 
Huid whirl vibrations do not exist, and the resulting self- 
excited vibrations are ofa whip type. This is not the case for 
the rotor/bearing system. The stiffness of the supporting 
springs, Ky. if ever installed, is usually small. The fluid film 
stiffness is even lower, so the last term in Eq. (6) can practi- 
cally be neglected in comparison to Кү. The rotor/bearing system 


instability threshold reduces, therefore, to 0 = т K | M 1. 
n r s Б 


This value may be much lower than the high eccentricity 
first balance resonance (as shown in Figure 2). А1 such a low 
speed, the post-instability threshold selFexcited vibrations 
arc fluid whirl, Note that the instability threshold £63 is 
inversely proportional to А. A lower A produces a higher 
instability onset, Here, therefore, is where one must look for 
the svstem stability improvements. Antiswirl injections [13. 
L4] arc an excellent way to increase stability. Note, also, that 
the onset of instability does not depend on the amount of 
Fluid damping, ГЭ), 

The remaining two cigenvalues from Eq.(3) provide two 
more natural frequencies of the system, The latter are clase 
toa the natural frequency of the forward and backward berd- 
ing modes of the rotor itself. Using Eq. (4), it can be shown 
chat, for the positive complex cigenvalue with the rotor nat- 


Gi (К, + К.) Мо, the imaginary рагі 


ural frequency 
reaches zero at a certain speed. This corresponds to the 
whip cessation ([17]; see Figure 62. 

The rotor limit сусе of sell-excited lateral vibrations of 
whirl ar whip, which occurs after the onset of instabilitv (6) 
can be found by solving Eqs. (3). The transition from a sta- 
ble rotor, r,2r;2ü into the limit cycle is not considered here. 
‚г;=А„ ci" (7) 

where A, and As are amplitudes of the rotor disk and jour- 
nal limit cycle circular precession, б), is the frequency of the 
whirl or whip.and б is the relative phase angle between the 


(b) 


rA, cor 


n Ct) 


r5 (t) 








Figure 5 
Hotor/seal (a) and rotor bearing (b) models, 
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disk and journal (or rotor at seal location? vibrations. These 
parameters can be calculated from two algebraic equations 
obtained by substituting Eqs 47) алго), After some transfor- 
mations, they are as follows: 


(K +K = Mey + Јо о = A£2) - 


| (8) 
+K,+K,+K,+K,(A.)|- Ky =0 
де = К,А, {K +K, -Mo + Јо (9) 


From Eq. (85, the Frequency, б), and amplitude, As, can be 
calculated. Eq. C9 provides A, and of, Ec. 68» has virtually the 
same format as the characteristic equation (4) with s=fop, 
except for an additional nonlinear term Kalaa). (Note that 
with the assumed general form of the Fluid film radial stiff- 
ness nonlinearity, this function now becomes the algebraic 
function of the amplitude, A.) The frequency of the self- 
excited vibrations, б), will be close to one of the natural fre- 
quencies calculated from Eq. (4). Therefore, it will be either 
close to A (for whirl or the rotor high cccentricity natural 
frequency of the bending mode, ,, (for whip). 

Eq. (©) provides a clue which allows identification of the 
axial location of the instability source of whirling or whip- 
ping in the machine train rotor. The phase, œ, is the relative 
phase between the whirl or whip vibration at the instability 
source and the other portion of the rotor. Since qo, the fre- 
quency of Fluid whirl or whip, is generally lower than wp: 
the angle € is, therefore, negative, Therefore, the disk sell- 
excited vibration lags the journal one, where the instability 
occurs. By measuring the phases of the filtered whirl or 
whip vibrations at two or more axial locations along the 
rotor the instability source can, therefore, be pinpointed, 
This method is now successfully used in the industry [1%]. 
When the rotor is at whip, the Frequency, e is close to Uh, 
and g= HOP, The rotor mode at whip is, therefore, a 
“corkscrew shape. 


Higher modes of fluid whirl and whip 

Extensive experiments, including modal perturbation 
testing and the identification of rotor system modal parame- 
ters, have produced more accurate models of Вот с] 
forces in seals and lighth-loaded bearings. These experi- 
ments have also led to the discovery and explanation of 
undocumented earlier phenomena, such as higher-smode 
fuid whirl (Figure б) and higher mode whip sel-excited 
vibrations [16-18]. 


Physical factors affecting fluid whirl and whip 

In this article, we wall discuss several physical factors that 
influence the fluidanduced, selbexcited vibrations (whirl 
and whip), 


Radial Sideload Force 

А rotor unbalance, generating the radial rotating force, 
may temporarily stabilize the rotor by suppressing the whirl 
during its гопар or coastdewn in the balance resonance 
speed: range. The application. of higher unbalance is not, 
however, an advised сиге to avoid Fuid-induced vibrations, 
constant, rotative speed independent, radial sidcload force 
applied to the rotor is more effective; This force can be gen- 
erated by rotor alignment, by process fluid distortion in 
Huid-handling machines, or finally, by gravity. when the jout- 
nal side of the horizontal rotor is made heavier! In the rotor 
rig (Figure L5, the radial force can be controlled by support- 
ing springs. The radial sideload forces made the journal 
move toa higher eccentricity inside the bearing, At the high- 
er eccentriciby, two positive factors operate: a larger fluid 
film stiffness and a lower A (see Figure 33. as the predomi- 
nantly circumferential flow pattern. changes, By increasing 
the direct bearing stiffness, the mechanical resonance fre- 
quency is increascd, By lowering A, the fluid-induced reso- 
nance frequency is decreased, Therefore, higher rotor 
eccentricity increases the margin of stability bv separating 
and Fluid-induced resonance, 
Both these factors affect the stability threshold (Oy An 
increase of Ky and а decrease of A result in higher £2, 
hopefully above the operating speed range. 


the mechanical. resonance 


А word of caution should be added here regarding the 
applied racial force magnitude: a small amount of the force 
may sufficienrlv stabilize the rotor ба “friendly siceload”), 
Too much farce would worsen bearing lubrication condi- 
tions and impose a stress on the rotor, which, ultimately, can 
leac ta cracks and the rotor breakage Gan “unfriendly side- 
load). 


Rotor Confipureatio 

Traditionally the cure Jor fluid-induced instabilities was 
sought in redesigning the rotor by making it shorter, fatter, 
and more rigid. This results in a higher first natural Freuen- 
cy, and in a higher instability threshold (6), as Kj and Ke 
increase, Analyzing, Eq. (б), vou will note, that for the rotor 
bearing system (Fig. 5b), the major stiffness component of 
the instability threshold is Kj. For a constant total stiffness 
K+ Ko and a constant 4,3 small К means that the rotor mass 
is very close to the journal. It also means that the instability 
threshold occurs at lower rotative speed.A high К, improves 
the stability situation. An immediate conclusion is to move 
the major mass of the rotor away from the fluid dubricatec 
bearing and maintain required journal eccentricity at the 
bearing using other means than gravity. This fact was con- 
firmed experimentally [3]. 
Bearing Geometry 

For a long time, bearing designers have known that a 
plain cvlindrical bearing is the worst design for rotor stabili- 
ty. Improvements were achieved by moditving bearing 


L Mme, however that making thee nor heavier dees nor mean tharian impyrevement is achieved. From once aspect, this is due te the existence or other radial fences, such is supe 
poring spring forces (if the springs exista сит journal ссу ае, Peon аги ис" йж it is eue pun tse noble of rotor mass ole im Er ig the higher nass ch is the fener 


instability theeshald is! 
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geometry. Elliptical, lemon-shaped bearings, bearings with 
grooves and lobes, and bearings with mobile parts, such as a 
tilting pad bearings, are generally more stable than plain 
cylindrical bearings. These improved features are associated 
with changes in the circumferential flow, which, in all these 
improved designs, has a fower strength, duc to boundary 
conditions, than in cylindrical bearings. Therefore, A is small- 
cr anid the instability threshold is (65 higher The decision to 
replace machine bearings with à morc stable design should 
be, however, justified by the analysis that the bearings (not 
the main flow) are the machine problems. 
Rotating Bearing Around a Stationary 
Post, Rough fournal Surface 

Interesting experimental results were obtained when 
the bearing. shell was mounted on the shaft, and the jour- 
nal was simulated by a stationary post. When maintaining 
the same geometrie dimensions and lubrication pattern as in 
the traditional rotor/bearing system (Figure 13, the whirl fre- 
quency was 0.5142 (as opposed to 0.470 for the regular 
bearing). Therefore, the increased linear velocity, A 
(Reradius) of the outer rotating surface, which drags the 
Fluid. inte rotational motion, increased the Fluid circumfer- 
ential average velocity [107]. In other experiments, a steel 
journal became rusty, duc to the use of an inappropriate 
lubricant. [ts acquired surface roughness and higher friction 
caused an increase of the Fluid circumferential average 
velocitv ratio. and, consequently, ап increase їп whirl fre- 
quency from 0.47 £2 to 0.4962. 
Fluid Pressure 

For partially lubricated bearings (180° bearings"), the 
lubricant pressure should be controlled. For example, it 
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should be lower than 20 psi (958 Pa), as a higher pressure 
mav change the lubrication pattern from partial to full. This 
will create a fully-developed circumferential Flow, and morc 
favorable conditions lor whirl and whip. In the case of full 
lubrication, increased. fuid pressure, which proportionally 
increases the fluid radial stiffness, is beneficial for stability. 
Externally-pressurized (^hwedrostatic"? bearings are well- 
known for their excellent stability features. In seals, a higher 
leakage flow may lower the circumferential Flow (A), while 
deteriorating machine efficiency, may prevent Muic-induced 
vibrations, 
Fluid femperatire 

Traditionally, the first quick-fix tried when a machine 
rotor started whirling or whipping was ro change bearing 
lubricant temperature. Oil temperature changes were саву 
to accomplish, being external to the online operating 
machine. The direction of temperature changes was, howev- 
er, more difficult te determine. In approximately 40% of 
cases, cooling the lubricant resulted in rotor stabilization. In 
approximately 3055 of cases, heating stabilized the rotor, In 
the remaining cases, this fix was ineffective. Heating or cool: 
ing the lubricant mainly affects its viscosity, decreasing or 
Increasing the fluid radial damping, /2. This equally. atfects 
the stabilizing agent( in Eqs. (33) and the instability driving 
tangential force, DAG. Therefore, che result cannot be pre- 
dicted, unless the instability threshold versus lubricant tent- 
perature is analyzed separately. This relationship is virtually 
nonlinear In some cases, an increase of temperature and 
lower oil viscosity may also affect fluid film stiffness, which 
moves the journal to higher eccentricity. 


Heating lubricant in only one bearing may bring some 
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Figure 6 


Spectrum cascade of rotor vertical vibrations measured by 
transducer B with timebase waveforms and “orbits” pro- 
duced by inboard and outboard transducers А and B at rota- 
tive speeds 2 Керт and 10.2 Керт. The “orbits” show that 


I4 Orbit 


during whirl of the first mode, the rotor vibrates in phase, 
behaving as a rigid body, while during the second mode 
whirl the rotor ends are about 90° out of phase [15]. 
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additonal positive effects; Due to thermal expansion, the 
rotor alignment state can be positively affected, providing a 
friendly sidelead. In general, however, vou shouldn't expect 
to bring about a long-lasting cure for whirl, and especially 
whip vibrations, by changing the lubricant temperature. 

Once morc caution regarding lubricant temperature 
changes may be useful, A higher temperature may increase 
bearing clearance, which makes the bearing more prone to 
instability. It mav alse affect the readings of the proximity 
transducer, if its scale Factor is sensitive to temperature 
changes. 
Antiswirl Injections 

Fluid antiswirl injections made tangentially to the seals or 
bearings in the direction opposite to rotation are by far most 
cfficient way to stabilize à rotor [13.14]. This technique may 
bc used in a passive or an active control mode (Figure 4) 
[14]. 


Formal derivation of fluid forces 

The formal derivation of Fluid dynamic forces was based 
on an assumption of a short cylindrical bearing or seal, pres- 
surized at the ends, lubricated by an incompressible fluid 
with fully developed, laminar circumlerential flaw [19], A 
new solution of the Reynolds Equation was obtained in the 
lorm of an infinite series. When calculating fluid forces, 
stresses (not pressures!) were integrated over the displaced 
rotor surface (not the bearing surface!) Fluid forces in the 
rotor bearing Cor seal) system were then obtained assuming 
rator rotation and precession at an eccentric position inside 
the bearing Cor seal, The equations allow lor extracting the 
elements of the Fluid force (13. 


Closing remarks 

This article summarizes the results of experiments on 
Fluid-nduced rotor sellexcited vibrations, whirl and whip. A 
simple mathematical model of the fluid force at low ro medi- 
um eccentricities, which was identified through modal per- 
turbation testing [4-6] adequately predicts the rotor 
dynamic behavior The rotor/bearing/seal model can also be 
extended to include other pertinent factors, such as fluid 
inertia. Fluid radial stiffness, damping anisotropy and damp- 


ing nonlinearity [2, 6], rotor lateral nonsvmmcetry, and more 
rotor modes. The key issue here is to adopt the multimode 
system (clastic solid body plus fluido approach in modeling. 
The rotor scli-excited whirl or whip vibrations can occur at 
апу mode, acquiring а frequency corresponding to the svs- 
tem natural frequencies. 
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